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We examine a Peierls ground state and its competing metastable state in tire one-dimensional quarter- 
filled Peierls-Hubbard model with the nearest-neighbor repulsive interaction V and the electron-phonon 
interaction (o= 1 /K with K being the elastic constant). From the mean-field approach, we obtain the phase 
diagram for the ground state on the plane of parameters V and K. The coexistent state of the spin-density 
wave and the charge ordering is realized for large V and K. With decreasing K, it exhibits a first-order phase 
transition to the unconventional Peierls state which is described by the bond-centered charge-density-wave 
state. In the large region of the Peierls ground state in the phase diagram, there exists the metastable state 
where the energy takes a local minimum with respect to the lattice distortion. On the basis of the present 
calculation, we discuss the photoinduced phase observed in the (EDO-TTF)2PF5 compound. 
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1. Introduction 

Electron correlation in quasi-one-dimensional molecular 
conductors has been studied extensively''^ where the spin 
density wave (SDW) state with the momentum 2kp (kp de- 
notes the Fermi momentum) originates from the combined 
effect of repulsive interaction and the nesting of the Fermi 
surface. The former effect is relatively large in molecular 
conductors and the latter effect becomes perfect for the one- 
dimensional band. In particular, the quarter-filled band sys- 
tems, which have been the main subject in molecular con- 
ductors, supply a rich variety of electronic states,^ e.g., the 
charge ordering (CO)'*'^ with a periodic array of charge dis- 
proportionation (h-1, 0, H-1, 0, • • • )■ In addition to the pure 2kp 
SDW, the coexistent state of the 2kp SDW and 2kp charge- 
density wave (CDW) has been observed in the X-ray exper- 
iment on (TMTSF)2PF6,^'^ where the coexistent state has a 
purely electronic origin due to the absence of lattice distor- 
tion. It has been clarified theoretically that the next-nearest 
repulsive interaction gives rise to such a coexistence,**'^ and 
that the state undergoes a first-order phase transition into the 
normal state with increasing temperature T.'" 

The molecular compound (EDO-TTF)2PF6, which is an- 
other salt showing a typical quasi-one-dimensional system, is 
a recent topic of the photoinduced phase transition. ' '"'^ It has 
been reported that the first-order phase transition from metal 
to the Peierls insulator occurs at T ~ 278 K and is followed 
by the fourfold periodicity of the charge-rich sites and the 
charge-poor sites with the spatial variation of (+1, 0, 0, +1, 
■ ■ ■). The large hysteresis is accompanied by the lattice distor- 
tion indicating the strong coupling to the lattice through the 
electron-phonon (e-p) interaction. This Peierls state shows a 
distinctive feature of the 2kp CDW where the amplitude of 
the lattice distortion takes a maximum between the hole-rich 
sites of the EDO-TTF molecules (i.e., the bond ordering). Al- 
though the spatial variation of the charge density is similar 
to the conventional CDW,^"'" the spatial pattern of the lat- 
tice distortion found in (EDO-TTF)2PF6 has not been fully 
understood in the context of the previous theory.'^ Further- 



more, the conductivity of (EDO-TTF)2PF6 exhibits a gigantic 
photoresponse, and the ultrafast insulator-to-metal transition 
is induced by the weak laser photoexcitation. The mechanism 
of this behavior has been discussed by assuming the existence 
of a metastable state without lattice distortion. However, it re- 
mains unclear if such a metastable state can be understood 
using the Peierls model in the presence of the electronic cor- 
relation. 

In the present paper, we examine the possible origin of the 
metastable state, which arises from the interplay of the Peierls 
state and the CO state. In §2, the model with both the repul- 
sive interactions and the e-p interaction is described and the 
formulation is given within the mean-field theory. In §3 and 
§4, the ground state and the metastable state are calculated 
to obtain the phase diagram. The possible parameter region 
for the metastable state is estimated. Section 5 is devoted to 
summary and discussions. 

2. Formulation 

We consider the 1/4-filled Peierls-Hubbard Hamiltonian 
given by'^' 

H = -tl^ £ (c]^c,+i.a + H.c.) 
j J 

+ f^M,(c]_^Cy+Ka + H.C.) 
jo 

j odd; 

where cj^a (spin a =T,i) is the annihilation operator of the 
electron at the site j and rij = Y.a a'-J-o- We impose the pe- 
riodic boundary condition, cj+^.a = Cj.a- Quantities U and V 
denote the on-site and nearest-neighbor-site repulsive interac- 
tions, respectively. The term proportional to uj (correspond- 
ing to the lattice distortion) denotes the e-p interaction and 
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the last two terms represent the elastic energy with the elastic 
constant K and K + 8K for even and odd j, respectively. The 
alternating elastic constant given by 5K plays a crucial role in 
the present paper As for the lattice distortion Uj, we take into 
account only the modulation with the momentum 2kf{= ^/2) 
(corresponding to the lattice tetramerization), which is rel- 
evant to the Peierls state in the (EDO-TTF)2PF6, and Uj is 
rewritten as 



MtCOS 



(2) 



where Ui is the amplitude of the lattice tetramerization. The 
elastic constant K is scaled so as to include the e-p cou- 
pling constant. In the present paper, we do not consider the 
possibility of the lattice dimerization with the momentum 
4A:f (= n)}^' Its detail is discussed in §5. 

After applying the Fourier transform, we define order pa- 



rameters (m = 0, 1,2,3) as''' 
1 



>mgo 



E E Sgn((7) {cl_„Ck+,nQa.a 



a=\,[-n<k<K 



MF 



_ 1 ^ 



'^k.a'^k+mQo.a 



MF 



(3) 



(4) 



c='\.[-K<k<n 

where Qq ~ 2kp ~ k/2, Sq = 0, Dq = 1 /2, and 
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We note that ^i, S2, and D2 (which are real numbers) 
correspond to the amplitudes for the 2kp SDW, 4kp SDW, 2kF 
CDW, and 4kp CDW, respectively, and 9 and 9' are the phases 
for the 2kp SDW and 2kp CDW. It has been found that the 
relation 9' = n/2 + 9 holds for — 0.^ In terms of these 
order parameters, the mean-field Hamiltonian is written as 



^^MF =E| (^* + ^+^) 4,a^k,a+ ^Ag^.^, + ^MtC'^ (1 - i)(cOS^ - sin^) j C^^g^^^Q^^ + H.C. +^2Qo,acl„Ck+2Qg,a^ 



N 
4 



U[-+2Di-2Si + Di-S'2]+Vil-4Di) 



(5) 



where e^^ ~ —2tcosk and 

U , M 

AQo,a = Y£>eo-sgn((7)-%, (6) 

Azeo.a = (I -2^) ^220 -sgn(cT)|52eo. (7) 

The total energy is given by E{ut,£,) = {Hmf)mf/N, where 
(• • • )mf is the expectation value over the mean-field Hamilto- 
nian. The ground-state energy is obtained by minimizing the 
total energy with respect to ut and ^ , where the corresponding 
equations are obtained as 

= e'^^-!^E(cos^-sinfc)(c;^Q+go^„^^^ + c.c., (8) 

k.o 

- — Mtsm(2^) 

= ^'^ ^ E('^°s^ - sin^) ^c|^Q.+e„,„^^^ + C.C.. (9) 

li.a 

The parameters Di, 02,81,82, ut, 9, 9', and ^ are evaluated 
self consistently using eqs. (3), (4), (8), and (9). We also ex- 
amine E{ui,^) to investigate the possible metastable state at 
Mt = 0. The present calculation has been performed for f/ = 3, 
where we take the hopping energy t as unity. 

In Fig. 1, we present the states obtained in our Hamiltonian. 
Two kinds of phases of ^ and 9 correspond to those of order 
parameters for the lattice tetramerization and the 2kp SDW, 
respectively: 

MtCOS (|; + <^), 5ie'^ 

The symbol * in Fig. 1 implies the state with Ux = 0. The SDW 
+ CO {D2+8\) is the charge-ordered state with the alternation 



of the charge-rich/poor sites. The pure SDW {81) has a spin 
amplitude, the maximum of which is located on the bonds in 
order to gain the transfer energy. Both the SDWh-CO state and 
the SDW state do not have any lattice distortion, namely, ui = 
0. The bond-centered CDW and site-centered CDW states are 
the Peierls states, i.e., the nonmagnetic states with finite lattice 
distortion (tetramerization) and with 2kp charge modulations. 
The lattice tetramerization takes the maximum amplitude on 
the bonds in the former state, while it does so on the sites in 
the latter state. We note that all these states exhibit the insu- 
lating behavior due to the presence of the 2kp density wave 
{81 or Di), which yields the gap in the dispersion at the Fermi 
energy. 

We note the relevance of the present analysis to the state of 
the quasi-one-dimensional (EDO-TTF)2PF6 compound. The 
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Fig. 1. Several possible ordered states obtained in the present model. The 
notation , 9) denotes the set of the phases for the lattice tetramerization 
{Uj = HtCos{;ry/2 + ^)) and for the order parameter of the SDW (5ie'*) or 
the CDW (Die'('"/-+'')). The ciirle with the shaded region denotes charge- 
rich sites and the bold line coiresponds to the bond with large transfer 
energy. 
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Peierls state observed at low temperatures is precisely of the 
bond-centered CDW type.^' Therefore, in the present paper, 
such a type of the CDW is mainly examined to find the possi- 
ble metastable state, which could be the origin of the photoin- 
duced phase realized after releasing the lattice distortion. 

3. Effect of Alternating Elastic Constant on Ground 
State 

In order to clarify the starting point of the present work, we 
begin to examine the case for 5K = 0, in which some results 
can be compared with the previous works. 

First, we briefly recall the electronic state in the absence of 
the e-p interaction.^ There is a critical value Vc ~ 0.28 where 
the state with 5i ^ and Z)2 = is obtained for V < Vc, and 
that wifli Si^O and Z)2 7^ is obtained for V > K. The CO 
state for V > Vc has a charge disproportionation with two-fold 
periodicity (Fig. 1). The phase transition at V = Vc is of the 
first order where the both order parameters D2 and Si exhibit 
discontinuous change. The difference in SDW between 
the case of V < Vc and that of V > Vc is that the maximum of 
the spin amplitude is on the bonds in the former case while it 
is on the sites in the latter case. Such a difference can be un- 
derstood from the commensurability energy for the Ikf SDW 
(Sq^ = 5ie'^) at quarter filling, which decreases continuously 
and vanishes only at V = Vc. Actually, the 6 dependence of 
eq. (5) takes the form'^ 

Eg{e) = const. +C{V)cos4e, (10) 

where the V dependence of C(V) is given by C(V) = Coo ^ 

Cqi V — CqsV^ -I with Coj being positive numbers. Thus, 

one finds C(V) > 0, i.e., 6 = n/4 (bond-centered density 
wave), for V < Vc, while one finds C(V) < 0, i.e., = (site- 
centered density wave), for V > Vc. We note that, due to such 
a sign change of C(V), the collective mode for the phase fluc- 
tuation exhibits the vanishing of the commensurability gap at 
V = Vc leading to the metalhc property. In fact, a noticeable 
behavior emerges in the optical conductivity g{(o) ((O being 
the frequency), where the static conductivity (7(0) becomes 
finite at V = Vc even for the cormnensurate SDW state.^*^ 

Next, we consider the case with the e-p coupling but with 
5K — 0. In Fig. 2, we show the phase diagram of the ground 
states obtained from eqs. (3), (4), (8), and (9). The bold 
Une denotes the boundary for the first-order phase transi- 
tion, which is estimated by comparing the minimum energy 
of E{ut,^). For large K (i.e., for small e-p coupling), we ob- 
tain Ut = Q and there is a critical value Vc ^ 0.28, where the 
SDW slate is obtained for V < Vc and the SDW+CO state is 
obtained for V > Vc. With decreasing K, the state with finite 
lattice tetramerization (Mt ^ 0) appears. There is the Z)i -|-5i 
state in between the SDW state and the bond-centered CDW 
state, where both boundaries do not depend on V. Such a V- 
independent result is ascribed to the mean-field treatment in 
which V does not contribute to the mean field for Si and Di 
[see eqs. (6) and (7)]. For the intermediate K and small V, 
there exists the region of the bond-centered CDW. This state 
resembles the BCDW state suggested by Clay et al. (Fig. 3(b) 
and Fig. 4 in ref. 18) and the DMh-SP state suggested by 
Kuwabara et al. (Fig. 2 in ref. 15), in the sense that the am- 
plitude of the lattice tetramerization takes a maximum on the 
bonds; i.e., bond-centered state. 



1.2 




0.5 1 1.5 2 



V 

Fig. 2. Phase diagram on the plane of V and K for 5K = and f7 = 3. 
The bold (thin) solid Une corresponds to the first- (second-)order phase 
transition. The region framed by the dotted rectangle is enlarged to show 
clearly several kinds of phases. 

Here, we note the effects of the quantum fluctuation and 
the lattice dimerization, which are not taken into account in 
the present mean-field treatment. It is possible that the mag- 
netic state of the SDWh-CO turns into a nonmagnetic state in 
the presence of the quantum fluctuation and the e-p couphng. 
Compared with the case in ref. 15, the mean field stabilizes 
the CO and the charge disproportionation in which the strong 
correlation and the quantum fluctuations take important roles. 
Thus, in the mean-field treatment, the boundary between the 
bond-centered CDW state and the site-centered CDW state is 
shifted to the weak couphng region, and the SDWh-CO state 
is stabilized for relatively small V. It is also expected that, 
due the quantum effect, the boundary between the SDWh-CO 
and site-centered CDW states becomes a crossover and the 
nonmagnetic state would be realized in both regions. Further- 
more the lattice dimerization can be expected to coexist with 
the bond-centered CDW state. '^''^ In the sense of the vari- 
ational principle, the region for the bond-centered CDW is 
extended by introducing the lattice dimerization. 

The bond-centered CDW state originates from the facts that 
the on-site repulsive interaction U separates two electrons 
being on the same site, and that the bond-centered ordering 
gives rise to the large gain of the transfer energy. For large 
V, the site-centered CDW state is reaUzed in order to avoid 
the increase in the energy of V. As seen from Fig. 1, the site- 
centered CDW state takes a maximum of the charge on the 
sites, where the gain of the transfer energy is made favorable 
owing to the electron hopping into both nearest sites. This 
state would correspond to the Akp CDW-SP state suggested 
by Clay et al. (Fig. 3(d) in ref. 18) and the CO+SP state sug- 
gested by Kuwabara et al.}^ since the lattice tetramerization 
takes the maximum amplitude on the sites. We note that, in 
the present case, the maximum of charge is located on the 
same position as that of lattice tetramerization (i.e., the am- 
plitude of Di is larger than that of D2). However, for the state 
in ref. 15, the maximum is obtained at the same position as 
that of the Ak-p CDW (i.e., the amplitude of D2 is larger than 
Di). This may originate from our choice of relatively small U 
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Fig. 3. Phase diagram on the plane of SK and K for U = 3 and V = I. The 
narrow region for small SK and around ~ 1 .03 is characterized by Z)i + 
Di + Si (see Fig. 2). With increasing SK, the region Di (bond-centered 
CDW) extends while that for Di +D2 (site-centered CDW) diminishes. For 
SK 7^ 0, the Di +D2 state takes ^ with an intermediate value of — w/4 < 
^<0. 



and then our site-centered CDW state is expected to move to 
the COh-SP state for large U. 

Now we show how the site-centered CDW (D\ +D2} di- 
minishes when 5K increases from zero. In Fig. 3, the varia- 
tions of the D\+D2 state and the bond-centered CDW (Di) 
are shown as a function of 5K. The SDWh-CO {D2 + Si) re- 
mains unchanged since Mt = 0. The addition of the 8K term 
increases the energy of the Di +D2 state with ^ ^0, but does 
not affect on the Di state with ^ = 0. Thus, the change from 
the D1+D2 state into the D 1 state occurs with increasing 5K. 
In the phase diagram of Fig. 4, the boundaries between the 
D\ +D2 state and the Di state are shown for the choices of 
5K = 0.05 and 0.1. The region of bond-centered CDW in- 
creases rapidly with increasing dK, i.e., a small amount of 
5K is enough to obtain the bond-centered CDW state as the 
ground state. It is noticed that the boundary of the direct tran- 
sition from the bond-centered CDW state to the SDW-nCO 
state emerges even for 5K = 0.05. Such a boundary is in con- 
trast to that in Fig. 2. 

Here, we compare the experimental findings in the (EDO- 
TTF)2PF6 compound with the results of the present calcu- 
lation. The Peierls state observed in (EDO-TTF)2PF6 is of 
the type of the bond-centered CDW,'^ and the photoinduced 
phase transition from the Peierls insulator to a metal takes 
place by the weak laser photoexcitation. This behavior has 
been discussed by assuming the metastable state without lat- 
tice distortion. ' ' In our calculation, the bond-centered CDW 
state is actually reproduced in Figs. 2 and 4. One can ex- 
pect the metastable state in the region near the boundary be- 
tween the the bond-centered CDW state and the undistorted 
state, since the the phase transitions are of the first order. 
Thus, we focus on the properties of the first-order phase tran- 
sitions from the bond-centered CDW state to the SDW-i-CO 
state. From the inset of Fig. 2, it can be seen that the system 
does not show a direct transition between the bond-centered 
CDW state and the SDW-nCO state, but there are interme- 
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Fig. 4. Phase diagram on the plane of V and K with several choices of 
SK = 0.05 and 0.1. The phase boundaries for finite SK are shown by the 

dotted line (SK = 0.05) and the dashed line (SK = 0.10). The region of the 
bond-centered CDW (Di) state is enlarged with increasing SK. For SK ^ 
0, the state Di +D2 takes ^ with an intermediate value of — w/4 < ^ < 0. 



diate states between these two states. Actually, the follow- 
ing three types of transitions occur when the bond-centered 
CDW state moves to the SDW-nCO state with increasing K. 
(i) For 0.67 < V < 1 .06, the bond-centered CDW state shows 
a first-order phase transition to the D1+D2 + Si state and un- 
dergoes a second-order transition to the SDWh-CO state, (ii) 
For 0.49 <V< 0.67, the bond-centered CDW state shows 
a second-order transition to the Di + Si state followed by the 
successive transitions of a first-order one into the Di -I-D2 -|- 5i 
state and a second-order one into the SDWh-CO state, (iii) 
For 0.28 <V < 0.49, the bond-centered CDW state shows 
a second-order transition to the Di + Si state and then un- 
dergoes a first-order transition to the SDW+CO state. The 
first-order phase transition in case (i) can be ascribed to the 
difference between the symmetry in the bond-centered CDW 
(E, =0) and that in the Di +D2 + S1 state (£, = -7t/4), where 
both states have finite lattice distortion. The origin of the first- 
order transition in case (ii) is also the same as that in case 
(i). In these cases, the second-order transition between the 
D1+D2 + S2 state with ^ = -7r/4 and the SDW + CO state 
impUes that the state around Mt = is unstable and turns out 
to be irrelevant to (EDO-TTF)2PF6. For case (iii), the first- 
order transition is due to the difference in the locking of the 
phase 6, i.e., 6 = n/4 in the Di +Si state and = in the 
SDW+CO state. This case is also irrelevant to the state of 
the (EDO-TTF)2PF6, since the spin ordering is absent in the 
Peierls state. Thus, both the metastable state at Mt = and 
the ground state of the bond-centered CDW (ui ^ 0) cannot 
be explained within the model of dK = 0, suggesting a new 
mechanism for the photoinduced phase observed in the (EDO- 
TTF)2PF6 compound. Our idea is the introduction of the alter- 
nation of the elastic constant, which enables us to reproduce 
the experimental findings. By increasing 5K, the region of the 
bond-centered CDW state is enhanced and moreover, there is 
a direct first-order phase transition to the lattice-undistorted 
SDW+CO state, as seen from Fig. 4. Such a behavior comes 
from the fact that the 5K has a role of fixing the phase ^ = 0, 
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which can be obtained from eq. (5). On the basis of such a 
consideration, we examine the bond-centered CDW state and 
the mechanism of the emergence of the metastable state at 
Mt = by setting ^ = in the next section. We discuss the 
relevance to the (EDO-TTF)2PF6 compound in §5. 

4. Phase Diagram and Metastable State for Large 5K 

In the following calculation, we examine the case of = 
corresponding to the bond-centered lattice distortion, which 
is reaUzed by choosing a moderate magnitude of 5K > 0.2 
as shown in the preceding section. The neglect of the lattice 
dimerization can be justified by considering the large 5K, as 
will be discussed in §5. 

In order to understand the role of the e-p interaction in the 
correlated system, we first examine eq. (5) by treating ux as an 
external field, i.e., by discarding the condition of eqs. (8) and 
(9). We calculate the order parameters Sj and Dj (j = 1, 2) as 
a function of Mt. The CDW, whose order parameter Di is al- 
ways finite due to the presence of the external field Ut, is called 
the extrinsic 2kp charge-density wave (ECDW). In Fig. 5, we 
explain how the state varies with increasing wt. When V =0, 
the order parameter Si decreases monotonically and becomes 
zero at Ut ~ 0.5. For V = 1, both the order parameters D2 
and Si take finite values for ut = 0, but as ut increases, D2 de- 
creases rapidly and becomes zero at ut ~ 0.26 and Si becomes 
zero at ut « 0.5. TheDi +Si state is understood by noting that 
Si is rather robust compared with D2 because U > V. With 
increasing ut (0.28 < V), the phase of the SDW, 6, increases 
from zero and reaches 7r/4 at which D2 vanishes. Such a vari- 
ation of 9 comes from the competition of D2 and Di. We ob- 
tain d' = d + n/2 for the Di + Si state and the Di state. For 
larger values of y(> 1.5), the strong competition of D2 and 
Mt results in the first-order transition, which is followed by the 
discontinuous changes of Sj and Dj {j = 1,2). Furthermore, 
for y (> 1 .7), the vanishing of D2 and Si takes place simulta- 
neously due to the strong effect of ut. 

The phase diagram is summarized in Fig. 5. The second- 
order and first-order transitions take place at the boundaries 
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Fig. 5. Phase diagram on the plane of ut and V where ut is regarded as an 
external field (see text). The thin and bold lines correspond to the bound- 
aries of the second-order and first-order transitions, respectively. Order pa- 
rameters Si,Di,D2 represent SDW, ECDW and CO, respectively. 



of the thin and bold lines, respectively, where the point with 
Mt = and V = Vc{— 0.28) is the singular one showing 
the first-order transition of D2. For the coexistent state of 
D\ +D2 + S1 +S2, both the amplitudes S\ and D2 are larger 
than that of Di. In the Di +Si state, the amplitude of Si is 
larger than that of Di . For large Mt, the pure ECDW (Di) state 
appears, e.g., for Mt > 0.5 and V < 1.7. The boundary between 
the Di -|- 5i state and the Di -|- £)2 + 5i + ^2 state in the region 
of 0.5 < y < 1.5 is given by y ~ 4Mt, indicating a direct com- 
petition between Mt and V. On the boundary corresponding 
to the vanishing of D2, one sees a point where the second- 
order transition changes into the first-order transition with in- 
creasing y. This variation resembles that found in the phase 
diagram in the V versus temperature (instead of Mt) plane.'" 
Such a fact may be explained by the existence of the term 
proportional to D2S^ in the energy expansion with respect to 
the order parameter'" The solid line of the phase boundary 
between the Di state and the Di + D2 + Si + S2 state is given 
by Mt 2± |y - 0.8 for y > 1.7. We found that Di is always 
finite due to the external field Mt, but is strongly suppressed 
in the presence of D2 (not shown), i.e., the strong competi- 
tion between D2 and Di, which is a characteristic of 1/4-filled 
systems. Also note that = ;r/4 if D2 = 0; i.e., the maximum 
density is located on the bonds in the absence of CO. 

Next we examine the Mt dependence of the mean-field en- 
ergy E{ut), which is defined by £'(Mt,<^ = 0). The minimum 
of this energy gives the true ground-state energy. In Fig. 6, 
the energy difference, 8E(ui) = E{ui) —£(0), is shown as a 
function of Mt for y = 0.1 with some choices of K. This be- 
havior denotes the conventional second-order phase transition 
except for its having several kinds of order parameters. Si and 
Di. For K = 1.06, the minimum is obtained at Mt = corre- 
sponding to the pure SDW state. The closed circle at Mt ~ 0.5 
denotes the point where the Si state vanishes. Thus, we ob- 
tain the SDW (5i) state for K > 1.05, the Di +51 state for 
1 03 < < 1 .05, and the pure bond-centered CDW {Di) state 
for K < 1.03. Figure 7 shows the case for y = 1. The novel 
feature in SE(ut) compared with that of Fig. 6 is the emer- 
gence of a local minimum at Mt = 0. For 0.79 < 1.01, we 
find a metastable state where the energy of the local minimum 
^(Mt = 0) is larger than the energy of the true minimum at fi- 
nite Mt. This fact implies that, with decreasing K, the system 
exhibits the first-order phase transition from the SDW-hCO 
(£)2 + Si) state into the bond-centered CDW state at « 1 .01 . 
We note that, for 0.28 < 0.70, there is a very narrow re- 
gion of K around K Ri 1 .03 for the Di -|- Si state (not shown), 
as found in Fig. 4. For a larger choice of y = 1.6, the local 
minimum at Mt = is found for 0.6 < 0.89. 

Here, we examine the origin of the metastable state found 
in Fig. 7. The mean-field energy of eq. (5) can be expanded in 
terms of Mt as 

E{ut)=E{0)+ (^a + ^K^u^ + Put + --- . (11) 

The coefficients, a and j3, can be estimated numerically and 
are shown in Fig. 8 as a function of V. With increasing V, 
there is a jump in both a and j3 at y = yc(~ 0.28), at which 
the pure SDW state moves to the SDWh-CO state; i.e., the D2 
state emerges. It can be found that a < for arbitrary V, and 
that j3 < (j3 > 0) for V >Vc(V < Vc). The result of a < 
is quite reasonable when we note that Mt acts as the external 
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Fig. 7. Mt dependence of SE(ut) for V = 1 with fixed K = 0.70, 0.95, and 
1 . 1 . With increasing Ht , D2 (Si ) decreases and vanishes at the point shown 
by the open (closed) circle. 



field. In order to understand why /3 is negative and diverges as 
y — > Vc + 0, we examine eq. (5) in terms of the phase degrees 
of freedom 0.'*' Actually, from eq. (5), the energy with fixed 
Mt and 9 is expressed as 

E{ut, 6) =Eg{d) + (a' - 7sin20 + ^K^ 



+ J3V + 



(12) 
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Fig. 8. V dependence of a and JS in eq. (11), where = ;r/4 (0) for V < Vc 
(V > Vc). The inset denotes — l/j5, which is proportional to (V — Vc) for 
V>V.. 



and 7(> 0) are estimated numerically from eq. (5). The J- 
term, which is proportional to 5j, denotes the interaction 
between the bond-centered CDW (Di) state and the SDW 
{Si) state. The coefficients are estimated as a' ~ —0.197, 
^' = 0.0068, and 7 ~ 0.056 for V = 1, where J3' ~ 0.0145 
for y < Vc. There is a small jump of 7 at V = Vc, where the V 
dependence of 7 is much smaller than that of a. The locking 
position of 9 can be determined in order to minimize eq. (12). 
The commensurability potential Eg{9) for V > Vc favors the 
locking position 0=0, while the 7-term favors 9 = 7t/4. For 
V > Vc, by minimizing eq. (12) with respect to 9, we obtain 
eq. (11) with the coefficients 



a = a 



I5=li'- 



(13) 



where £^(0) is given by eq. (10). Quantities a'{< 0), /3'(> 0), 



s\c{v)\- 

We can immediately reproduce the anomalous behavior 
-1/J3 oc (y - y^) atV = Vc+ 0, by noting that the coeffi- 
cient C(y) follows C(y) tx: (y - Vc). This is the reason why 
the coefficient of the fourth-order term j3 can become negative 
for y > yc. It is shown that the result similar to eq. (13) can 
be obtained from the bosonization scheme, in which the spin 
degree of freedom is also taken into account.'^' For V < Vc, 
the 7-term does not contribute to the wf term since both the 
commensurability potential and the 7 term favors the lock- 
ing position 9 = n/4 (one simply obtains a = a' — f). Thus, 
the conventional second-order transition is reproduced. From 
these arguments, the conditions for the metastable state are 
summarized as follows: (/) the existence of the CO for y > yc 
(the commensurability energy leading to 9 = 0), (ii) the cou- 
pling between the SDW state (Si) and the bond-centered CDW 
state ( Ut), which is given by the 7 term in eq. ( 12), and ( Hi) the 
mechanism to fix ^ =Qas discussed in §5. 

Finally we show the phase diagram of the ground state on 
the plane of V and K in Fig. 9. For K > 1.05 correspond- 
ing to the weak e-p coupling, we obtain either the pure SDW 
(y < Vc) with ;?r/4 or the SDWh-CO (V > Vc) with = 0. 
With decreasing K (i.e., increasing the e-p couphng), the 
bond-centered CDW state with 9 = n/4 appears for K < Kc 
where Kc is a critical value of the first-order phase transition 
for the Peierls state. The state with ut = Q becomes metastable 
in the large interval region ofK<Kc as shown by the shaded 
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Fig. 9. (Color online) Phase diagram of the bond-centered CDW (Di), 
SDW+CO (D2 + 5i ), SDW (5i ), and bond-centered-CDW+SDW (Di + 5i ) 
states, on the plane of V and K. The shaded region denotes the area where 
the metastable state at Wt = appears. The thin and bold lines correspond 
to phase boundaries of the second-order and first-order transitions, respec- 
tively. 

area. The lower boundary of the shaded region is given by 

K=-Aa, (14) 

which is verified from a in Fig. 8. There is also another 
metastable state with the local minimum at Ux ^ 0, which is 
located in a certain region of K > Kc (not shown in Fig. 9). 

5. Summary and Discussion 

We examined the Peierls-Hubbard model at quarter-filling 
with the intersite interaction (V), and obtained the ground- 
state phase diagram within the mean-field theory. A no- 
ticeable finding is the undistorted state iui = 0), which be- 
comes metastable in the bond-centered CDW state close to 
the boundary between the bond-centered CDW state and the 
SDWh-CO state. The obtained bond-centered CDW ground 
state and its metastable state may share the following common 
features with the Peierls state of the (EDO-TTF)2PF6 com- 
pound. The spatial variation of the Peierls state is the same, 
and the metastable state for Ui — Q could be related to the pho- 
toinduced phase found in the experiment on the EDO-TTF 
compound. 

In the present analysis of §4, we assumed ^ —Q in order 
to stabilize the bond-centered CDW state and examined the 
metastable state at ui = 0, which may be relevant to the prop- 
erty of the EDO-TTF compound. Actually, in the normal state 
of this compound, there is the alternation of the bending of the 
molecule for every two sites along the the one-dimensional 
chain. In addition to the ground state (§3), we discuss this 
assumption for the metastable state on the basis of the energy, 
which is expanded in terms of Ui. The energy with fixed Mt, Q, 
and can be exphcitly written as 

E{ut,0,^) =C(y)cos40 + (a'-7sin(2e-2<^))M^ 

+ (f + ^sin2^)«? + ..., (15) 



where the 5K term denotes an increase in elastic energy for 
1^ 7^ 0. Here, we note that, even for V > Vc, there is no 
competition between 9 and ^ for dK = 0, since the mini- 
mum of eq. (15) is obtained by choosing ^ so as to satisfy 
sin(20 -2^) = 1. In this case, the y^/C{V) term in j3 [see 
eq. (13)] is absent and the energy ^(Mt, 9,^) [eq. (15)] takes a 
true minimum at ut = (i.e., j3 > 0). When 5K ^ 0, the energy 
E{ut,9,£,) as a function of £, increases and the metastable 
state around ut = Q can be expected (i.e., j3 < 0). Actually, for 
dK > 0.46, the energy E{ut, 9,£,) takes a local minimum at 
Ut — for all values of in the case of K — 0.9 and V — 1. 
Thus, the assumption ^ = can be justified when the degree 
of alternation of the elastic constant 8K becomes larger than 
a critical value. 

In the present analysis, we did not consider the possibility 
of lattice dimerization (uj = (— l)-'Md, where is the am- 
plitude of the lattice dimerization). Since the elastic energy 
of the lattice dimerization is given by {N/2){K + 5K/2)u^, 
5K also has an effect of suppressing the lattice dimerization. 
From these arguments, we expect that the effect of 5K will 
arise from the bending freedom of the molecules in the (EDO- 
TTF)2PF6 compound. Such a bending plays important roles 
in the photoinduced cooperative transition, although the ori- 
gin of the molecular property still remains an open question 
from a microscopic view point. 

Here, we note states at finite temperatures. There are some 
theoretical works on finite-temperature properties. Quite re- 
cently, it has been pointed out that the Peierls-Hubbard model 
with 8K = exhibits the first-order transition from the dimer- 
Mott state into the spin-Peierls state,^^ where both states show 
the insulating behavior On the other hand, the purely elec- 
tronic model shows the transition from the metallic state into 
the insulating state of the 2kp SDW + 4kp SDW + 2kp CDW 
state. Compared with the state of the (EDO-TTF)2PF6 com- 
pound, the normal state at high temperatures is different from 
the dimer-Mott state obtained from the former model and the 
lattice deformation cannot be reproduced in the latter model. 
It is expected that the parameters in the shaded region in Fig. 9 
will lead the first-order phase transition into the normal state 
with increasing temperature, and that such a transition will 
be strongly enhanced by the effect of the large bending of 
molecules in the (EDO-TTF)2PF6. 
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